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ABSTRACT 

The  viscoelastic  properties  of  dilute  aqueous  solutions  of 
polyethylene  oxide  were  studied  using  a  modified  simple  elastometer. 
The  dynamic  viscosity  and  dynamic  rigidity  are  calculated  from  the 
measurements  in  the  frequency  range  from  1  to  11  Hertz  for  concentra- 
tions of  .25%,  .5%,  1%,  and  1.5%  of  a  polymer  having  a  molecular 
weight  of  about  4  million. 

The  results  indicate  that  the  solutions  may  be  represented  as 
a  simple  Maxwell ian  element  in  the  frequency  range  considered. 
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1.      Introduction 

Polyethylene  oxide,    a  high  molecular  weight,    long-chain  linear 

molecule  with  the  structure   shown  below, ** 

Ay  \     y 

has   been  the   subject  of  many  recent    investigations.      When  dissolved   in 
water,    it   forms   a  viscoelastic   solution  which  exhibits   some  of   the 
characteristics  of  both  a  solid  and   a   liquid.      That    is,    it  has   an 
elastic  modulus   and  a  coefficient  of  viscosity.      Polyethylene  oxide 
solutions   also  exhibit   a  decrease   in  the  apparent  viscosity  with   in- 
creasing rates  of   shear,    which   is   a  characteristic  of   so  called  pseudo- 
plastic  or  non-Newtonian  fluids.      Philippoff   first   found  this   phenomenon 
for   certain  fluids;    he  also   found  that   the  viscosity  of  polymers   under 
oscillating  deformations   depends  only  on  the  frequency  and   not   on  the 
amplitude  of  deformation.    (7) 

Non-Newtonian  fluids   are  best   characterized  by  a  complex  vis- 
cosity or   impedance  function    il  —   «C  ""  •  'W.      Where  \C|    is  the   usual 
dynamic  viscosity  and    Vi—  J2,    where    (3     is  the  dynamic   rigidity  and 
^0    is   the  angular  frequency  of  oscillation. 

The  real   and   imaginary  parts  of   »L    give  good   insight    into  the 
properties  of   the   solution  under  periodic  deformation.      If  Y(j    is 
zero,   the  fluid   is  perfectly  viscous,    and   if    iC^  is    zero,    the  fluid 
is   perfectly  elastic.     The  use  of  a  complex  impedance,    however,    is 
useful  only  for  oscillatory  motion. 

Polyethylene  oxide  solutions  have  been  the   subject  of  two   thesis 
studies   at   this    institution.      Lt.    Chester  found,    in   1964,   that  these 
solutions   behave  essentially   like  water   for  sound  waves   at   relatively 


high  frequencies    (10  to   26  megaHertz) .    (1)      The   long-chain  molecules   in 
solution  contribute  no  measureable  effects   to  the  absorption  or   speed 
of   sound.      He  did   find,   however,    that   shear  waves    in  the  megaHertz 
range  are   strongly  affected  by  the   addition  of   the  polymer.      In   1965, 
Lt.    Kinnier  and  Lt.    Reister   found,    using  a  torsional  oscillator,    that 
from   118  Hertz   to   675  Hertz,   viscosity  changes   with   frequency,    and  that 
the   long-chain  molecules  have  a  great  effect  on  the  properties  of  the 
solution.    (4) 

The  object  of  this    investigation  was   to   determine  the  viscoelastic 
properties  of   the  solutions   in  the    low  oscillatory  frequency  range  from 
1   to    10  Hertz.      The  polyethylene  oxide  used    in  this   investigation  was 
purchased   from  the  Union  Carbide  Company.      It   is   named   Polyox  WSR-301; 
its   average  molecular  weight    is  approximately  4  X   10      and  the  distri- 
bution  is   unknown. 


2.   Equipment 

In  order  to  measure  the  complex  viscosity  at  low  frequency  in 
Polyox,  a  modified  simple  elastometer  was  constructed,  following  the 
designs  of  Goldberg  and  Sandvik,  (3)  and  Markovitz,  Yavorsky,  Harper, 
Zapas  and  DeWitt.  (6)  Goldberg  and  Sandvik  designed  and  constructed  a 
simple  elastometer  which  provided  a  quick  and  simple  method  for  deter- 
mining viscoelastic  properties  of  gels  and  soap  solutions.   Its  opera- 
tion depends  on  the  motion  of  a  cylindrical  bob  suspended  in  fluid 
undergoing  periodic  angular  oscillation.   Its  operation,  however,  is 
limited  to  the  frequency  of  natural  resonance  for  the  suspended  cyl- 
inder.  Markovitz,  et  aL  took  the  basic  design  of  Goldberg  and  Sandvik 
and  extended  its  frequency  range. 

The  equipment  used  in  this  investigation  is  similar  to  the  elasto- 
meter of  Markovitz,  but  utilizes  an  optical  method  for  obtaining  the 
needed  measurements.   It  consists  essentially  of  a  circular  cylinder 
containing  a  liquid,  in  this  case  Polyox,  in  which  another  smaller 
cylinder  is  suspended  concentrically.   The  outer  cylinder  is  made  to 
oscillate  mechanically  through  a  small  angle  with  a  frequency  that  can 
be  varied.   The  motion  of  the  inner  cylinder  with  respect  to  the  outer 
cylinder  can  be  used  to  determine  the  complex  viscosity  of  the  fluid. 

The  apparatus  (Figures  1  and  2)  consists  of  a  steel  table  with 
three  leveling  legs,  on  which  a  turntable  with  a  lathe-turned  shaft  is 
mounted.   The  turntable  is  supported  by  ball  bearings  and  is  rigidly 
connected  to  a  light  aluminum  rocker  arm.   A  smooth  rubber-covered 
wheel  mounted  at  the  end  of  the  rocker  arm  is  held  by  a  spring  against 
an  eccentric  cam.   The  eccentric  cam  is  driven  by  a  shaft  and  pulley 
system  coupled  to  a  variable  speed  electric  motor.   Since  the  drive 


motor  is  more  stable  in  operation  at  high  speeds,  two  sets  of  pulleys 
are  used  to  reduce  the  shaft  speed.   The  pulley  system  also  isolates 
the  steel  table  from  any  motor  vibration. 

If  the  eccentricity  of  the  cam  is  small  compared  to  the  length  of 
the  rocker  arm,  then  the  motion  of  the  turntable  can  be  considered 
sinusoidal . 

A  metal  can  is  mounted  on  the  turntable  with  a  smooth  glass  cyl- 
inder cut  from  a  length  of  large  diameter  glass  tubing  fixed  in  the 
center.   This  glass  tube  serves  as  the  outer  cylinder  and  the  metal  can 
as  a  container  for  a  constant  temperature  bath. 

Several  different  diameter  inner  cylinders  were  constructed  using 
machined  and  polished  Lucite  rod.   The  inner  cylinder  is  suspended  by 
a  thin  (No.  38)  copper  wire  fixed  to  a  screwjack  device  which  permits 
the  suspension  height  to  be  varied.   The  screwjack  device  can  also  be 
moved  laterally  to  permit  centering  of  the  inner  cylinder. 

Optical  means  are  employed  to  measure  the  motion  of  the  inner  and 
outer  cylinders.   This  system  functions  as  follows:  a  mirror  is  mounted 
on  the  side  of  the  outer  cylinder  and  another  mirror  is  mounted  on  the 
inner  cylinder.   These  two  mirrors  are  illuminated  with  a  highly  colli- 
mated  beam  of  light  from  a  common  light  source.   The  reflected  beams 
of  light  are  displayed  on  a  curved  screen  placed  three  meters  from  the 
inner  cylinder.   The  curved  screen  eliminates  any  distortion  of  the 
end  points  of  the  arcs  of  reflected  light.  Although  the  reflecting 
mirrors  are  not  the  same  distance  from  the  screen,  the  distance  of  sep- 
aration is  small  compared  to  the  distance  to  the  screen  so  that  the 
error  caused  by  this  is  negligible. 

Two  photocells  are  mounted  behind  slits  on  the  curved  screen. 
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They  can  be  positioned  in  the  midpoint  of  the  arc  of  reflected  light 
from  the  oscillating  cylinders.   The  output  of  each  one  of  these  cells 
is  amplified  and  displayed  on  separate  traces  of  a  dual  beam  cathode 
ray  oscilloscope  (Figure  3).  This  permits  measurement  of  the  phase 
angle  between  the  displacements  of  the  inner  and  outer  cylinder  by 
using  a  vernier  delay  dial  to  align  the  output  traces.   The  amplitude 
of  the  oscillation  is  measured  directly  on  the  face  of  the  screen  by 
measuring  the  length  of  the  light  arc. 

There  was  a  slight  amplitude  modulation  of  the  magnitude  of  both 
reflected  arcs.   This  may  have  been  caused  by  a  minor  wobble  in  the 
turntable,  resulting  from  the  light  load  on  the  supporting  thrust 
bearings. 

This  apparatus  is  useful  in  the  frequency  range  from  1  to  11 
Hertz.  At  higher  frequencies  the  rocker  arm  tends  to  float  or  vibrate 
on  the  eccentric  cam,  and  at  lower  frequencies,  the  size  of  the  inner 
cylinder,  the  restoring  force  of  the  wire,  and  their  relation  to  the 
viscosity  of  the  fluid  becomes  very  critical. 
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3.      Theory  of  Measurement 

The  measuring  device   (Fig.    4)    consists  of   a   lucite  cylinder  of 
radius  r   suspended  to   depth  h   in  a   liquid  contained   in  another  cylinder 
of  radius  R.      The  outer  cylinder    is   caused  to   execute  simple  harmonic 
oscillations   through  a  small   fixed   angular  amplitude.      The  relationship 
between  the  motion  of  the   inner  cylinder  and  the  viscoelastic   proper- 
ties of  the  liquid  was   developed  by  Markovitz.    (5,6) 

Consider  a  cylinder   lamina  of  the   liquid  of  height  h  with   inner 
radius   z  and  outer  radius   z   +  dz.      The  torque  which  acts  on  the  cylin- 
drical  surface  of  the  cylindrical   lamina  at   z   is    —  l.T\  ™     |Ch^(*®L 
where   <&    is   the  angular  velocity  of   the   liquid  at    z  and  \t    is  the  vis- 
cosity coefficient,    which  may  be  either  real   or  complex.      At   z   +  dz, 
the  torque  on  the  surface   is  ^H^Y?  If)  g®\  *h  \Ot&fMM  ]  J-  »    and  the 
net   torque  on  the  lamina   is  ^ 


£  [in  e*rih  &*)]}* 


As   a  result,   the  element  of  volume  will   acquire   an  angular   acceleration 
Q     such  that^ng  h/>  ©  J^  Z  &  fett  £*h  YUlH/  3  '**£    Since  the  motion   is   a 
forced  oscillation  with  frequency  ^~-     the  substitution      (3     -  {p\  £ 


can  be  made,    and  the  equation  of  motion  becomes 

The  solution  of  this  equation   is        Qi    |  i)  Ji  fo£)  lp  1  ii^*)J    where  J 

— » 
and  Y  are   Bessel    functions  of   the  first  order,   A  and  B  are  arbitrary 

'     h 
constants   to  be  determined  by  the  boundary  conditions   and  C^    is    (-J (jji'_ \ 

At  the   inside  wall  of  the  cylinder  containing   the   liquid, 

JOJC 
the  motion  can  be  represented   by       ©($~^  Oft  ^       since  the   liquid   is 
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following  the  forced  motion  of  the  cylinder  which  has  amplitude  (S  f% 
the  suspended  lucite  cylinder  experiences  a  shearing  torque  equal  to 
£FT  V*  fm  f*®  1  due  to  the  liquid  and  a  restoring  torque  —  K6H- due  to 


ta^l 


the  suspension  wire  which  has  a  torsional  constant  K.   As  a  result  the 
bob  acquires  an  angular  acceleration  such  that 


where  I  is  the  moment  of  inertia  of  the  cylinder.   Putting  the  boundary 
conditions  into  the  solution  for     one  can  arrive  at  the  equation  of 
motion  for  the  inner  cylinder. 


r^r      = 


gib  3tyr3  [J~X^)X(^  "IM  X^)J 


(*  ^)[U^)M~W^^ 


where  fc*  - jUU_h  JQ  W       QnC(   J^^j   ^  f^V^e  0/  dx(0?)  ^^r  k 

Although  this  represents  the  exact  relationship  between  the  desired 
quant  ity\D  and  the  measured  quantities   ^   it  is  in  an  inconvenient 
form  since  ((  ,  the  unknown,  is  buried  in  the  complex  argument  of  the 
Bessel  functions.   Making  a  Taylor  expansion  of  the  functions  with 
argument  C4.  V-  about  Q[^  the  following  result  is  obtained: 
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6y 


+  j_ 


I^-K        R    -Jlf     +    cu/°    (R^-H) 


a  L  Wfth  co  R>* 


8 


ft 


a* 


(rI^)(/?w-srV.2^)+iirV^j?\ 


nno 


/h  - 


Ai  - 


6>  = 


I 


^TThR'r 


Jw-J. 


2L    U^+  4't'I 


c^  * 


r 


and  neglecting  higher  order  terms  results  in 


ev- 


er 


-,[coS/,^j,  ,  fi  (fc^-i]  -ij(ft„^ 


Yc  now  may  ^e  found  by  the  solution  of  a  complex  quadratic  equation  or, 

J. 
if  the  term  involving  ,-  ^may  be  neglected,  then  by  a  first  order  equa- 


tion. 


] 


h***^  ~X  $M?    *X   faifftMo*  -6 


iyVi 


ZJ 


rO 


In  this   investigation   it  was   possible  to  neglect  all   terms  of 

I 

higher  order  than  "T  ,  because  of  the  errors  in  the  measurement  of 

It 
phase  angle  and  amplitude. 

All  of  the  constants  may  be  calculated  from  direct  measurements 

except  for  the  torsional  constant  of  the  suspension  wire.   The  torsional 

constant  may  be  determined  from  the  frequency  of  oscillation  of  the 

inner  cylinder  suspended  in  air  using  the  relation 

The  values  of  the  constants  for  the  three  inner  cylinders  used 
are  shown  in  Table  1. 
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4.  Experimental  Procedures 

Basically,  the  operation  consisted  of  filling  the  inner  cylinder 
with  the  sample  of  Polyox  to  be  tested,  then  immersing  and  centering 
the  inner  cylinder.   The  relative  magnitude  of  oscillation  and  phase 
lag  angle  is  then  recorded  as  a  function  of  oscillation  frequency. 

The  samples  of  Polyox  were  all  mixed  in  a  2000  ML  beaker  with  a 
three  bladed  beater  at  100  RPM.   The  powdered  Polyox  was  slowly  added 
during  the  mixing  cycle  and  care  was  taken  not  to  form  clumps  in  the 
solution. 

The  turntable  was  carefully  leveled  before  a  run  to  insure  that 
the  inner  cylinder  would  be  suspended  concentrically  with  the  outer 
cylinder. 

The  inner  cylinder  was  centered  before  the  Polyox  was  introduced 
by  checking  its  position  over  a  pointed  plug  which  fitted  in  the  outer 
cylinder  center.   The  inner  cylinder  then  was  lowered  to  a  prescribed 
depth  after  the  outer  cylinder  had  been  filled. 

The  response  of  the  elastometer  is  very  dependent  on  the  viscosity 
of  the  fluid,  the  distance  of  separation  between  the  cylinders,  and  the 
frequency  of  oscillation.   As  an  example,  in  a  1%  solution  of  Polyox 
with  a  12  mm.  diameter  inner  cylinder  and  a  14.4  mm.  diameter  outer 
cylinder  at  frequencies  of  less  than  3  Hertz,  both  cylinders  appear 
to  move  together,  and  at  frequencies  higher  than  8  Hertz,  the  motion 
of  the  inner  cylinder  becomes  very  small.   Accurate  measurements  of 
viscosity  are  difficult  under  either  of  these  conditions.   The  radius 
of  the  inner  cylinder,  or  the  restoring  force,  must  be  varied  in  order 
to  cover  the  frequency  range  of  interest.   In  this  investigation,  dif- 
ferent inner  cylinders  were  used  with  the  same  outer  cylinder. 
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Most  of  the  problems  encountered  were  mechanical  in  nature,  such 
as  the  centering  of  the  inner  cylinder,  and  the  stability  of  motor 
operation. 

Although  the  apparatus  was  constructed  with  a  constant  tempera- 
ture jacket,  this  was  not  employed  during  the  runs.   The  measurements 
were  made  at  room  temperature,  which  was  24  _  1°C.   The  dynamic  vis- 
cosity of  Polyox  is  temperature  dependent  but  Chester's  results  indi- 
cate that  any  uncertainty  caused  by  a  *   1°  temperature  variation  would 
be  less  than  the  experimental  error. 
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5.   Results  and  Discussion 

The  results  are  best  shown  graphically  (Figures  4,5)  with  the 
plots  of  K|  the  dynamic  viscosity,  and  sr  the  modulus  of  rigidity, 
vs.  frequency  for  different  concentrations  of  Polyox.   Selected  values 
are  also  presented  in  Table  2.   These  two  quantities  are  plotted  with 
a  straight  line  approximation  for  different  concentrations.   The  dynamic 
viscosity  decreases  with  increasing  frequency  of  oscillation  while  the 
modulus  of  rigidity  increases. 

\c,   and  Cz      are  also  plotted  (Figure  6)  as  a  function  of  concen- 
tration at  a  frequency  of  3  Hertz.   Both  f[,   and  C3  change  greatly  with 
increasing  concentration  and  both  of  the  plots  have  the  same  general 
shape. 

An  exact  estimation  of  the  amount  of  error  is  difficult  to  make  in 
this  investigation.   The  elastometer  was  tested  with  a  Newtonian  solu- 
tion of  known  viscosity  and  the  results  checked  within  10%.   The  error 
appeared  to  be  random,  so  it  could  not  be  used  as  a  corrective  factor 
to  be  applied  to  other  measurements. 

The  most  critical  variable  was  the  phase  lag  angle.   Since  \L\ 

and  (j    are  direct  functions  in  the  computation  formula  of  the  sine 

and  cosine  of  the  phase  lag  angle,  they  were  markedly  affected  by  any 

errors  in  that  angle.   This  was  particularly  true  when  the  angle 

measured  was  near  O.  2L  or  Tf   where  the  sine  or  cosine  change  most 

2 
rapidly. 

The  phase  angle  is  determined  by  measuring  the  ratio  of  the  sweep 

delay  to  the  period.   For  small  delays  the  uncertainty  in  measurement 

due  to  lack  of  trace  stability  and  eye  estimation  is  of  the  same  order 

of  magnitude  as  the  observed  quantity.   An  estimate  of  the  uncertainty 

in  the  phase  angle  is  1  3°. 
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A  phase  angle  meter  was  also  used  to  try  to  improve  the  accuracy 
of  the  measurement;  however,  it  was  not  capable  of  giving  good  results 
on  the  irregular  pulse  outputs  from  the  photocells. 

The  amplitude  radio  was  easily  measured  to  within  an  accuracy  of 
t  5%  since  both  light  beams  swung  through  an  arc  ranging  from  1%  to  8 
inches  in  length. 

In  a  given  series  of  measurements  with  a  particular  inner  cylin- 
der, the  amplitude  ratio  starts  at  unity  for  low  frequencies,  increases 
to  a  maximum  at  system  resonance,  then  decreases  as  frequency  increases, 

The  phase  lag  angle  starts  at  zero,  increases  to  ^    at  resonance,  and 

2 

approaches    as  the  frequency  increases.   A  value  of  \ii   and  <j      was 
determined  at  whatever  frequency  the  amplitude  ratio  and  phase  lag 
angle  was  measured,  and  these  values  were  used  to  plot  the  dependence 
on  frequency.  An  alternate  method  would  be  to  plot  the  phase  lag 
angle  and  amplitude  ratio  as  a  function  of  frequency  and  use  values 
determined  from  this  plot  to  calculate  [C,  and^  .   Regardless  of  the 
method  used,  each  point  on  the  plot  would  have  a  different  estimate  of 
error  due  to  uncertainties  in  determination  of  frequency,  of  phase  lag, 
of  amplitude  ratio,  and  of  the  other  physical  constants  of  the  system. 
The  points  with  the  largest  estimated  error  are  the  end  points  in  a 
particular  series  of  measurements,  i.e.,  those  at  the  low  and  high 
frequencies.   The  phase  lag  angle  for  these  points  is  near  O  and  JJ" 
and  any  uncertainty  in  it  affects  the  results  greatly. 

The  straight  lines  were  plotted  by  judgment  and  an  overall  error 
of  _  10%  is  considered  reasonable,  based  on  the  results  of  measure- 
ments of  fluids  of  known  viscosity. 

While  the  investigation  was  in  progress,  it  was  learned  that 
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Michael  C.  Williams  at  the  University  of  Wisconsin  had  made  similar 
measurements  on  1.46%  Polyox  solution  using  a  Birnboim  concentric 
cylindrical  pumping  viscometer.   His  measurements  covered  a  wider 
frequency  band  but  his  results  match  the  results  of  this  investigation 
and  those  of  Lt.  Kinner  and  Lt.  Reister  very  closely.   His  results  are 
plotted  with  the  results  of  this  investigation. 

The  simplest  way  of  accounting  for  the  mechanical  properties  of 
Polyox  solutions  is  to  introduce  a  single  viscosity  iCo  and  a  single 
shear  modulus  (y&  .   One  method  of  superimposing  these  properties  is 
the  Maxwell  model.   In  the  Maxwell  model,  the  fluid  is  represented  by 
a  Maxwellian  element  consisting  of  a  spring  (representing  the  shear 
modulus)  and  a  dashpot  (representing  the  viscosity)  in  series.   If  \l$ 
is  the  zero  shear  rate  viscosity,  ^  is  the  frequency  of  the  applied 
disturbance,  then  the  Maxwell  model  for  viscoelastic  solutions  gives 
a  formula  for  the  dynamic  viscosity 

where 

T-    1  IT  ft  6  . 

(y<.&    -   the  shear  modulus  at  infinite  shear. 

Using  values  of  iCi  determined  in  this  investigation  at  3  Hertz 
and  6  Hertz  to  solve  the  equation  of  the  simple  Maxwell  model  for  ylo 
and  7   for  various  concentrations  of  Polyox,  values  for  '£&  are  com- 
puted which  are  the  same  order  of  magnitude  as  the  low  shear  rate 
viscosity  coefficients  measured  at  6  RPM  with  a  Brookfield  viscometer. 
The  values  are  compared  with  the  relaxation  times  in  Table  3.   The 
agreement  between  \(.0  is  better  at  the  lower  concentrations  and  these 
results  indicate  that  in  the  frequency  range  from  1  to  11  Hertz  the 
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behavior  of  Polyox  solutions  may  be  approximated  by  a  single  Maxwell 
element.   The  results  of  Lt.  Kinnier  and  Lt.  Reister  indicate  that  as 
frequency  increases,  this  approximation  is  no  longer  valid. 

The  elastometer  did  not  prove  to  be  a  good  experimental  device 
for  investigating  Polyox  solutions.   The  properties  of  Polyox  solutions 
change  greatly  with  frequency  and  this  requires  many  combinations  of 
inner  and  outer  cylinders  to  cover  a  limited  frequency  range.   This 
investigation,  I  feel,  has  filled  a  gap  in  the  Polyox  spectrum.   Fur- 
ther work  is  needed  in  smaller  concentrations  and  at  higher  frequencies. 
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Figure  3.  Block  Diagram  of  Measuring  Apparatus 
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